Optical surface modes in the presence of nonlinearity and disorder 
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We investigate numerically the effect of the competition of disorder, nonlinearity, and boundaries 
on the Anderson localization of light waves in finite-size, one-dimensional waveguide arrays. Using 
the discrete Anderson - nonlinear Schrodinger equation, the propagation of the mode amplitudes 
up to some finite distance is monitored. The analysis is based on the calculated localization length 
and the participation number, two standard measures for the statistical description of Anderson 
localization. For relatively weak disorder and nonlinearity, a higher disorder strength is required to 
achieve the same degree of localization at the edge than in the interior of the array, in agreement 
with recent experimental observations in the linear regime. However, for relatively strong disorder 
and/or nonlinearity, this behavior is reversed and it is now easier to localize an excitation at the 
edge than in the interior. 

PACS numbers: 42.25.Dd, 42.65.Wi, 42.79.Gn, 72.15.Rn, 73.20.Fz 



Introduction.- A fundamental question concerning sys- 
tems which are both disordered and nonlinear is whether 
or not Anderson localization jl| is weakened by the pres- 
ence of nonlinearity. While it was originally developed in 
order to understand electronic transport in non-periodic 
(disordered) soUds, the concept of Anderson localization 
was later generalized to the localization of classical waves 
in disorder media The interaction of propagating 
waves, and in particular of electromagnetic waves, when 
both disorder and nonlinearity are present can signifi- 
cantly affect localization and other phenomena Q. 

Despite of many efforts, that question has not been 
conclusively answered pl-[ll|. It thus seems that the an- 
swer depends on the relative strength of disorder and 
nonlinearity. For large nonlinearity, time-periodic and 
exponentially localized excitations in the form of discrete 
breathers may be generated, due to the self-trapping 
effect For small disorder strength, the discrete 

breathers are modulated to become localized modes [l^ . 
The above theoretical results were accompanied by a se- 
ries of experimental demonstrations of Anderson local- 
ization in optics [1^ and Bose- Einstein condensates . 

It was recently observed experimentally that Ander- 
son localization in finite segments of disordered waveg- 
uide arrays in the linear regime is actually site-dependent 
. Specifically, a higher disorder strength is required to 
achieve the same degree of localization at the edge than 
in the interior (i.e., the "bulk") of the array |16|]. Here 
we are interested in the effect of the interplay of disorder 
and nonlinearity on the site-dependence of wavepacket 
localization in one-dimensional (ID), disordered, finite- 
size arrays of coupled Kerr- type waveguides. Using the 
discrete nonlinear Schrodinger (DNLS) equation with di- 
agonal (on-site) disorder, frequently referred to as the dis- 
crete Anderson - nonlinear Schrodinger (DANLS) equa- 
tion, we calculate standard measures of Anderson local- 



ization in order to analyze the site-dependence of the de- 
gree of localization for a wide range of nonlinearity and 
disorder strengths. 

Model equations and statistical measures.- Consider a 
ID array of N single-mode optical waveguides. In the 
framework of the coupled-modes formalism, the electric 
field C{x,z) propagating along the waveguides can be 
expanded as a superposition of the waveguide modes, 
C{x,z) = ~ '^)' where C„ is the complex 

amplitude of the single guide mode (j){x) centered at the 
?ith site. The evolution equations for the modal ampli- 
tudes Cn are 

dC 

i—, \-Vn,n-lCn-l+^nCn+Vn.n+lCn+l+X\(^n\ Cn = 0, 

dz 

(1) 

where n — 1,2,3, N , e„ is the propagation constant as- 
sociated with the nth site, Vn,n±i are the tunneling rates 
between two adjacent sites, x is the nonlinearity parame- 
ter, and z is the spatial coordinate along the propagation 
direction ('time'). Eq. ([T]) describes very well recent ex- 
periments in ID disordered waveguide lattices and, more- 
over, it serves as a paradigmatic model for a wide class of 
physical problems where both disorder and nonlinearity 
arc important. Disorder is introduced into the optical 
lattice by randomly choosing the propagation constants 
e„ from a uniform, zero-mean distribution in the interval 
[—A, -I- A]. As a result, the lattice remains periodic on 
average and, to a very good approximation, the param- 
eters Vn.n±i become independent of the site number n, 
i.e., Vn,n±i — V. Then, Eq. ([T]) reads 

dC 

i—^ + e„C„ + V{Cn-i + Cn+i) + x\Cn?C^ = 0. (2) 
dz 

In order to take into account the termination of the struc- 
ture, we impose free boundary conditions at the edges. 
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i.e., Co = Cn+1 = 0. For x — ^ 0; Eq.© reduces to the 
original Anderson model while in the absence of disor- 
der (e„ — 0), it reduces to the ID DNLS equation 
that is generally non-integrable and it conserves the norm 
TV = J2n=i l^nl^ ^^'^ Hamiltonian "H. 

To investigate the simultaneous interplay of disorder, 
nonlinearity and boundary effects, we place initially a 
single-site excitation C„ = 5n,no near, or at the bound- 
ary of the array. This determines the value of the norm 
= 1 for all subsequent 'times'. For a quantitative 
analysis we utilize two of the standard measures used 
in the description of Anderson localization; the partic- 
ipation number P = ^J2n=i ' ^^'^ local- 
ization length £, defined as the width of the envelope 
containing the localized profile. The participation num- 
ber gives a rough estimate of the number of sites where 
the wavepacket has significant amplitudes, and it is a 
useful measure for ascertaining localization effects in the 
case of partial localization. In this case, P will saturate 
at a finite value, indicating the formation of a localized 
wavepacket. 

Statistical analysis.- In the following, we set V = 1, 
while the nonlinearity parameter x varies between and 
10, and the disorder width A takes on several different 
values. Since Anderson localization is essentially a sta- 
tistical phenomenon, many realizations of disorder are 
needed to obtain meaningful averages for the quantities 
of interest. This is particularly true for low-dimensional 
systems. We typically use hr = 1000 realizations in each 
run. The array contains N = 200 waveguides, and the 
maximum evolution "time" is z = 100 (except otherwise 
stated). In optics, nonlinearity is varied by changing the 
power content of the input beam. However, this is for- 
mally equivalent to keeping the norm of the wavepacket 
fixed, and to varying the nonlinearity parameter x- Eqs. 
([2]) are integrated with a standard 4rth order Runge- 
Kutta algorithm with fixed time-stepping. We compute 
the absolute squared profiles < |C„p >, where the brack- 
ets denote averaging over all realizations ur, hereafter re- 
ferred to as Anderson mode profiles. Assuming that the 
Anderson modes have a z— dependence which is a simple 
exponential function of the form < |C„p >= C^^^ e~^^^, 
the localization length £ can be computed via fitting 
procedure, with C^^^. being the numerically obtained 
maximum of < |C„p >. 

In Figs. 1-3 the inverse localization length is shown 
as a function of x for three different values of disor- 
der strength. In all the cases displayed in these figures 
the initial wavepacket is a single-site excitation placed at 
n = no, with no = 1 (right at the edge), 2, 3, 5, and 
10. In Fig. 1 (where A — 0.6) we easily identify two 
different x~ regimes; the weak and the strong nonlin- 
earity regime, where £~^ is small and large, respectively. 
That generaly implies a lower degree of Anderson local- 
ization in the weak nonlinearity regime compared to that 
in the strong nonlinearity regime. The large in the 
interval of x values where all the curves fall the one onto 




FIG. 1. (color online) Inverse localization length as a 
function of the nonlinearity strength x- A single-site excita- 
tion is launched from no = 1 (solid - black), 2 (red - dot- 
ted), 3 (green - dashed), 5 (blue - long-dashed), 10 (indigo - 
dotted-dashed), for a system with = 200, V ^ 1, z ^ 100, 
ur = 1000, and A = 0.6. 




FIG. 2. (color online) Inverse localization length and 
localization length £ (inset) as a function of the nonlinearity 
strength x- Same parameters and notation as in Fig. 1, but 
A = 1.0. 



the other, indicates the existence of a highly localized 
mode due to the self-trapping effect. The characteristic 
nonlinearity strength, Xcj that roughly distinguishes be- 
tween the two regimes is the critical on for self-trapping 
to occur in the ID DNLS equation fl^. In the weak 
nonlinearity regime another important feature appears; 
as it can be seen in the figure, the £~^{x) curve obtained 
for excitations initially placed at the edge (no = 1) is 
well below all the others (for which no > 1). Thus, an 
excitation initially placed at the edge (no = 1) leads to 
final wavepackets that are less localized than those which 
have been initialized below the 'surface' (no > 1). This 
effect can be understood as the "repulsive" action of the 
boundary, reported in a previous work for surface modes 
in nonlinear periodic lattices 0, and it is in agreement 
with the experimental observations of Ref. [l^ . 

As the disorder strength A is increased from 0.6 to 1.0 
(Fig. 2), all the curves become flatter without showing 
any qualitative difference from those of Fig. 1. When 
A is increased to 2.0, however, we do observe qualitative 
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FIG. 3. (color online) Inverse localization length and 
localization length £ (inset) as a function of the nonlinearity 
strength x- Same parameters and notation as in Fig. 1, but 
A = 2.0. 



differences (Fig. 3). For weak nonlinearity (x ~ 4 ^ 
Xc) there are no significant differences in the degree of 
localization for the Anderson modes resulting from initial 
excitations either at the edge or in the interior of the 
array. Thus, it is as easy to localize a wavepacket at the 
edge as it is in the interior in this case. However, for 

intermediate nonlinearities (from x ^ 4 to x ~ 8) it is 
more favorable to loacalize a wavepacket at the edge than 
in the bulk, whereas for large nonlinearities x ^ 8 it is as 
easy to localize a wavepacket at the edge as it is in the 
bulk. Thus, for relatively strong disorder we observe a 
bahavior that is strikingly different to what is observed 
in Figs. 1 and 2. The two different behaviors can be 
seen even more clearly by comparison of the localization 
length ^ as a function of x for A = 1 and 2, shown 
in the insets of Fig. 2 and 3, repsectively. Thus, the 
presence of strong disorder is capable of overcoming the 
" repulsive" character of the boundary for any value of x 
and, moreover, it favors wavepacket localization at the 
edges for intermediate nonlinearities. 

Typical examples of localized mode profiles both at, 
or close to the edge and the 'bulk' are shown in Fig. 4, 
where the Anderson mode profiles < jCnp > are shown 
as a function of n for A = 0.6. (Note that in this figure 
z = 1000.) The edge-localized modes are significantly 
more extended than the bulk modes, even though the 
former are not all localized exactly at the edge. This is 
because of the small disorder strength A, which allows 
the repulsive force of the boundary on the mode to dom- 
inate and move slightly the mode-maximum towards the 
bulk. However, similar profiles (not shown) are obtained 
also for A = 1.0, that is large enough to keep the local- 
ized modes at their initial location. 

Moreover, single-site excitations initialized at different 
sites no can be pushed by the boundary towards the in- 
terior and form Anderson modes at the same final site. 
These modes are different, at least for finite propaga- 
tion distances z; they differ in the degree of localization, 
leading to a multiplicity of Anderson modes having their 
maximum at the same site of the lattice (Fig. 5). For 




FIG. 4. Anderson mode profiles < |Cnp > as a function of 
the site-number n for = 200, 1/ = 1, A = 0.6, z = 1000, 
riR = 1000, and X = 1 (a,d), X = 2 (b,e) and x = 3 (c,f). 
Left panels denote the surface mode case (no — 1) while right 
panels refer to the bulk mode case (no = 100). Only part of 
the array sites are shown for clarity. 
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FIG. 5. (color online) Averaged absolute squared profiles < 
jCnp > as a function of the site- number n for N = 200, 
1/ = 1, A = 0.4, ur = 1000, X = 2.5, z = 100. These profiles 
result from a single-site initial excitation at no = 1 (black- 
solid); no = 2 (red-dotted); no = 3 (green-dashed). Only 
part of the array sites are shown for clarity. 



the particular value of x used for Fig. 5, three single-site 
excitations initialized at different tiq have formed, after 
they have been propagated up to z = 100, three distinct 
Anderson localized modes whose maximum is located at 
the same lattice site (i.e., at n = 3). 

Finally, let us look at the participation number P as 
a function of z for wavepackets that are initially local- 
ized at the edge (no = 1) and in the 'bulk' (no = 100) 
of the waveguide array. The logarithmic plots are shown 
in Fig. 6 and Fig. 7, respectively, for several values 
of X and weak disorder. Comparing the curves in these 
figures corresponding to the same x^ see that those 
for no = 100 are shifted to higher P values than those 
for no = 1. Thus, single-site excitations initialized at 
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FIG. 6. (color online) Logarithm of the averaged participation 
number logj^Q(< P >) as a function of logjQ(z) for the surface 
case (no = 1) and = 200, V = 1, A = 0.6, ur = 1000, 
z = 1000. The values of x are shown on the figure. 




FIG. 7. (color online) Logarithm of the averaged participation 
number logj^Q(< P >) as a function of logj^Q(2:) for the bulk 
case (no = 100) and = 200, y = 1, A = 0.6, ur = 1000, 
z = 1000, no = 100. The values of x are shown on the figure. 



no = 100 have, while propagating along z, a larger num- 
ber of sites where the wavepacket has significant ampli- 
tudes. However, the excitations initialized at tiq = 100 
exhibit a higher degree of localization than those initial- 
ized at no = 1 (see also Fig. 4). It is also interesting to 
see how the curves in each figure change as a function of 



X- For X = 6, well above Xc we see that the wavepacket 
remains localized at a single site, independently of ng. 
For X = 4 - Xc, the logio(< P >) vs. logio(2:) curve 
increases slowly with increasing log]^Q(z), and it satu- 
rates at a finite value around z ~ 100, indicating the 
formation of a wavepacket highly localized around uq. 
For the nonlinearity strengths that are less than Xc the 
logjg(< P >) vs. logjg(2:) curves exhibit qualitatively 
the same behavior. There is an increase with increas- 
ing log]^o(z) which is slowed down after some z specific 
to each x value, and indicates significant delocalization 
of the initially single-site wavepacket. Delocalization is 
stronger for decreasing nonlinearity strength. However, 
those curves do not seem to saturate, implying that the 
corresponding Anderson localized modes may delocalize 
further at longer propagation distances. 

Concluding remarks.- We have performed extensive 
calculations with the ID DANLS equation in order to 
clarify some aspects of the interplay between boundary 
effects, disorder and nonlinearity, in finite-size waveguide 
arrays. In particular, we attempt to clarify the site- 
dependence of Anderson localization that results from 
that interplay. We computed two standard measures of 
localization for discrete systems for varying nonlinearity 
and disorder strengths, and we observed two strikingly 
different behaviors depending on the strength of the dis- 
order. For weak to moderate disorder, we distinguish 
two different nonlinearity regimes; weak and strong, for 
values of x roughly below and above Xc, respectively. 
In the weak nonlinearity regime it is easier to localize a 
wavepacket in the interior of the array than at the edge, 
which is in agreement with the experiments in the lin- 
ear regime [l6|. In the strong nonlinearity regime it is 
as easy to localize a wavepacket at the edge as it is in 
the interior. However, for relatively strong disorder, this 
behavior is reversed, at least for intermediate nonlinear- 
ities, and it is now easier to localize a wavepacket at the 
edge than in the bulk. For weak and very strong non- 
linearities there is no significant site-dependence on the 
degree of wavepacket localization. The results presented 
here obviously hold for finite propagation distance z, an 
important case of practical interest for experimentalists. 
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